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Abstract 



Aims. To investigate the effect of gravitational lensing of supernovae in large ongoing surveys. 

Methods. We simulate the effect of gravitational lensing magnification on individual supernovae using observational data input from 
two large supernova surveys. To estimate the magnification due to matter in the foreground, we simulate galaxy catalogs and compute 
the magnification along individual lines of sight using the multiple lens plane algorithm. The dark matter haloes of the galaxies are 
modelled as gravitational lenses using singular isothermal sphere or Navarro-Frenk- White profiles. Scaling laws between luminosity 
and mass, provided by Faber-Jackson and Tully-Fisher relations, are used to estimate the masses of the haloes. 
Results. While our simulations show that the SDSSII supernova survey is marginally affected by gravitational lensing, we find that 
the effect will be measurable in the SNLS survey that probes higher redshifts. Our simulations show that the probability to measure 
a significant (3<r) correlation between the Hubble diagram residuals and the calculated lensing magnification is J; 95% in the SNLS 
data. Moreover, with this data it should be possible to constrain the normalisation of the masses of the lensing galaxy haloes at the la 
and 2cr confidence level with ~ 30% and ~ 60% accuracy, respectively. 

Key words, supernovae: general - gravitational lensing 



1. Introduction 

Type la supernovae (hereafter SNe la) are exceptionally use- 
ful tools for cosmological investigatio ns. Ongoing large su- 
pernova surveys, such as ESS ENCE ([M iknaitis e t alj l2007t 
IWood- Vasev etafll2q07l) . SNLS dAstier et al.ll2006h and SDSSII 
dFrieman et alj I2008I) harvest hundreds of distant supernovae 
with high precision. 

Several effects will alter the luminosity of these distant 
sources, and thus dilute the cosmological signal. In this pa- 
per we investigate gravitational lensing (de)magnification of 
the light from these explosions. Although this effect is hardly 
large enough to dominate the uncertainties in current exper- 
iments, as has already been addressed by several groups on 
statistical grounds (e.gjRiess et alj|2004t iHolz & Lindedl2"005r. 



the improved statistics may well compensate for the smaller dis- 
tances. To find out which surveys are most likely to display a 
lensing signal we have performed simulations based on our pre- 
vious work. The results are presented below, and show that the 
prospects of detecting lensing are very good. 

In Sect. [2] we briefly introduce the surveys we have inves- 
tigated and the data from these surveys that are used as input 
for the simulations. A short summary of gravitational lensing of 
SNe la is given in Sect. [3] Section [4] describes the simulations 
and the results are presented in Sect. |5J where we also discuss 
how a detection of a lensing signal could be used to obtain in- 
formation about the lensing matter. Finally, we summarise our 
results in Sect. [6] 



lAstier et al.l2006tlWood- Vasev et al j2007l) . the question if it can 
be corrected for on an individual supernova b asis remains. This 
has been investigated in a series of papers dGunnarsson et alj 
2006; iJonsson et al.ll2006h and a tentative detection of a correla- 
tion between calculated lensing magnification and Hubble resid- 
uals was recen tly found using the v ery high-z supernovae in the 
GOODS field dJonsson et alj2007t) . 

In this paper we investigate to what ex tent other surveys 
can be affected. The GOODS supe rnovae dRiess et all l2004t 
Strolg er et alj2004tlRiess et alj2007l) remain unique in that their 
large distances clearly make them more likely to be lensed. On 
the other hand, the ongoing ground based surveys will measure 
significantly more supernovae with much better sampling, and 
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2. Supernova surveys 

There are (at least) three ongoing major ground based supernova 
surveys as mentioned above. We have chosen to simulate the 
SNLS survey to estimate the effects in a deep survey. It probes 
a similar redshift range as the ESSENCE survey and is likely to 
find more supernovae. It is also of importance that SNLS uses 
several filters which allows photo-z determinations for the field 
galaxies. The SDSSII supernova survey is selected to investi- 
gate if any gravitational lensing effects could be expected for 
this smaller redshift domain when we have good statistics. 
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2.1. SNLS 

The Supernova Legacy Survey, (SNLS. IXstier et al.|[2Q06b has 
published constraints on the cosmological parameters using 71 
SNe la ranging from z — 0.2 — 1.0 discovered during the 
first year. Around 500 spectroscopically identified SNe la are 
expected to end up in the final Hubble diagram at the end of 
the survey. The survey consists of an imaging survey detecting 
and monitoring the light curves of the SNe, and a spectroscopic 
follow-up confirming the nature of the SNe and measuring their 
redshifts. A wide field imager used in "rolling search" mode al- 
lows to perform the discovery and photometric follow-up at the 
same time. The SNLS is searching for SNe in 4 different fields 
of 1 square degree each and each field is observed every third 
to fourth night for as long as it remains visible. Observations 
are taken in 4 filters (the MegaCam filter set) giving rise to a 
promising data set for observing the lensing signal. 

Since the line of sight to SNe la located close to the edge 
of the fields can not be modelled for gravitational lensing accu- 
rately, we estimate that only ~ 450 of the ~ 500 Hubble diagram 
SNe la can be used to search for a gravitational lensing signal. 
The parameters f rom this survey that are n eeded as input for this 
investigation (see Gunnarsson et al. 2006) have been taken from 
the literature (lllbert et al.ll2006l) and from discussions with the 
involved astronomers. 

2.2. SDSSII 

The Sloan Digital S ky Survey II Supernova Survey (SDSSIISN, 
Frie man et al.ll2.Q08h is searching a 280 square degree field over 
3 seasons. The survey is well on track to discover over 500 Type 
la supernovae, of which ~ 300 will likely end up in the Hubble 
diagram. The SDSSII is a follow-up on the previous success- 
full SDSS projects, which means that this search has a working 
infrastructure capable of handling large amounts of data. The 
search is performed in 5 bands (ugriz) and the SPSS fie ld is 
calibrated down to an accuracy of ~ 1% (llvezic et~al . 2003). 

3. Gravitational lensing of SNe la 

Here we give a short introduction to the aspects of gravitational 
lensing of SNe la which are relevant to this work. For a more 
t horou gh introduction to the subject see e.g. Bergstrom et al.l 
d2000h . 

Matter inhomogeneities in the Universe give rise to gravi- 
tational fields which will influence the path of light rays prop- 
agating between a source and an observer. In the weak lensing 
regime, gravitational bending of light, or gravitational lensing, 
can lead to changes in the apparent position of the source as 
well as changes in the flux. For point sources such as SNe la, 
the amplification (or de-amplification) of the flux can be de- 
scribed by the magnification factor, /i. The observed flux from 
a point source in an inhomogeneous universe is given by / Q b s = 
fif, where / is the flux that would be observed in a homo- 
geneous universe in the absence of lensing. Amplification and 
de-amplification relative to a homogeneous universe is conse- 
quently described by /j, > 1 and \i < 1, respectively. 

The distribution of magnification factors, for sources 

at a specific redshift depends on the distribution of matter be- 
tween the source and the observer. In general is asymmet- 
ric with a peak at fi < 1 and a high magnification tail. Due to 
flux conservation, as long as we do not have multiple images, 
the average magnification is unity, i.e. (//) = J P(/i)^G?/i = 1. 
The average flux from a number of standard candles, all at the 



same redshift, (/ bs) = (m)/ = /> i s trius an unbiased distance 
indicator if the variance of the flux is finite. However, changes 
in the flux from standard candles will increase the dispersion in 
distance measurements and could potentially lead to selection 
effects. 

4. Simulations 

We use simulations to investigate the effects of gravitational 
lensing for the surveys described in Sect. [2] Since gravitational 
lensing depends on the distribution of matter between the source 
and the observer, a model of the matter distribution is needed in 
order to simulate gravitational magnification. In the next section 
we briefly describe our model of the matter distribution in the 
Universe, which is based o n observational input. Fo r more de- 
tails we refer the reader to iGunnarsson et al.l (120061) where our 
model is described in detail. 

4. 1. Simulating lines of sight 

For each simulated SN la we simulate a line of sight populated 
by foreground galaxies. The galaxy popula tion is characterise d 
by the luminosity functions presented in Dahl en et al.l (|2005). 
These luminosity functions are used to simulate the number of 
galaxies along a line of sight as well as their brightness and type. 
Simulated B-band absolute magnitudes of the galaxies are in 
the range —23 < Mb < —16. A constant comoving number 
density of galaxies with redshift is assumed. Since the position 
of the galaxies in the sky as well as the redshifts are assigned at 
random, clustering of matter into larger structures is not taken 
into account. 

Each galaxy is associated with a dark matter halo and we 
assume that the mass of the halo can be obtained from the lu- 
minosity of the galaxy. Gravitational lensing depends not only 
on the mass of the lens, but also on the distribution of mat- 
ter within the lens. We will model the galaxy haloes using sin- 
gular isothermal sph ere (SIS) and Navarro-Frenk- White (NFW, 
iNavarro et al.| [l997) profiles. Since these profiles are divergent, 
they are truncated at r2oo, the radius inside which the mean den- 
sity is 200 times t he critical density of the Universe. We use nu- 
merical routines (INavarro et al.lll997h to find the concentration 
of a halo characterised by the mass m2oo enclosed within the 
radius r 2 oo- 

To convert galaxy luminosity to halo mass we relate dynam- 
ical properties of the galaxies to the dark matter halo. For a SIS 
halo, TO200 is related to the velocity dispersion, a, via the for- 
mula 

771200 = 5Glf ■ (1) 
This formula is a good approximation also for a NFW halo with 
the same corresponding maximum rotation velocity, V max . The 
mass of a NFW halo is over-estimated by 5 — 15% when Eq. (Q]) 
is used to calculate m 2 oo for a given value of V max . The veloc- 
ity dispersion of early type galaxies can be obtained from the 
absolute magnitudes through the Faber-Jackson (F-J) relation as 
derived bv lMitchell et al.l d2005l) 

logo- = -0.091(M B - 4.74 + 0.85z). (2) 

The redshift dependence in the relation accounts for the bright- 
ening of the stellar population with redshift. We represent the 
measurement error in this relation by the sc atter in the SDSS 
measurements reported bv lSheth et al.l d2003l) 

rms(logcr) =0.079[1 + 0.17(M B + 19.705 + 0.85^)]. (3) 
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For late ty pe galaxies we use the Tully-Fisher (T-F) relation de- 
rived by Pierc e & Tullyld 19921) . with correction for redshift cal- 
culated by Boh m et alj d2004l) . 

logK„a X = -0.134(M B + 3.61 + 1.22*0, ( 4 ) 

to relate galaxy luminosity to the maximum rotation velocity of 
the galaxy. The velocity dispersion relates to the maximum ro- 
tation velocity via a = V max / For the T-F relation the ob- 
served scatter is rms(Afs) = 0.41, which corresponds to 

rmsOogKiax) =0.06. (5) 

The mass range for early and late type galaxies computed using 
the formulae above is 4 x 10 11 < to 2 oo/-M© < 3 x 10 13 and 
3 x 10 10 < m 2 oo/M@ < 2 x 10 13 , respectively. 

The code used to compute the magnification factor for the 
simulated lines of sight is described in the following subsection. 

4.2. Q-LET 

The Quick Lensing Estimation Tool (O-LET. [Gunnarssonl2 004) 
is a numerical code which utilises the multiple lens plane algo- 
rithm. Each dark matter halo is projected into a lens plane per- 
pendicular to the line of sight situated at the angular diameter 
distance corresponding to the redshift of the galaxy. The deflec- 
tion angle is then computed for each lens plane and a ray orig- 
inating at the position of the image is traced back to the source 
position via the lens equation. From the Jacobian determinant of 
the lens equation, the magnification factor can be computed. 

In our model universe, the matter density of dark matter 
galaxy haloes, p g (z), is typically less than the measured global 
matter density p m (z) = £!mPc(1 + z ) 3 ' where p® is the present 
critical density of the Universe. Since only galaxies brighter than 
the magnitude limit will be associated with haloes, we expect 
p g (z) to decrease with redshift. To ensure that our model is con- 
sistent, we describe the matter which is "missing" due to our 
ignorance as a smoothly distributed component characterised by 
the smoothness paramete r given by rj (z) = 1 — p g (z) / p m (z0 

Numerical routines in lKavser et al.l(l 19971) . which can handle 
a redshift dependent smoothness parameter, are used to calculate 
the angular diameter distances in a clumpy universe which are 
involved in the gravitational lensing calculation^. 

4.3. Estimating magnification factors 

It should be possible to detect gravitational lensing of SNe la 
through the expected correlation between / D b s and p. However, 
since p is not observable we have to use an estimate, p cst , when 
we search for the correlation. To simulate the precision to which 
magnification factors can be estimated, we first compute p for 
a fiducial model and then try to recover the value of the mag- 
nification factor in the presence of different uncertainties. These 
uncertainties, which are described in the following paragraphs, 
dilutes the lensing signal carried by the estimated magnification 
factor p cat . 

Finite field Size Only a limited number of foreground galaxies 
can be taken into account in the calculations, since the fields are 

1 A completely smooth universe and its opposite, a completely 
clumpy universe, corresponds to T] = 1 and r\ = 0, respectively. 

2 The cosmological model used in the calculations is described by 
S1 M = 0.3, f2 A = 0.7, and H Q = 70 km s" 1 Mpc" 1 . 



of finite size. To limit this source of error to the percent level, 
all ga laxies within 60" should be included (Gunnarsson et alj 
2006). To include this source of error in our simulations a larger 
number of galaxies is included in the calculation of p (all galax- 
ies within 100") than in the calculation of p es t (all galaxies 
within 60"). 

Photometric redshifts Only galaxies brighter than the magni- 
tude limit can be included in our calculations of the estimated 
magnification factors, p cst . For surveys like SNLS we have to 
rely upon photometric redshifts in our calculations. The mag- 
nitude limit of the SNLS is i' AB = 25.5, but reliable photo- 
metric redshifts can only be obtained to i' AB = 24, which we 
take as the magnitude limit in our simulations. Photometric red- 
shifts, z p , usually have larger uncertainties than spectroscopic 
redshifts, z s . More accurate photometric redshifts can be ob- 
tained for bright g alaxies than for faint galaxies. According to 
Ilber t et aTl d2006l) the photometric redshift accuracy of galax- 
ies observed within the CFHT legacy survey with brightness 
17.5 < i' AB < 22.5 and 22.5 < i' AB < 24 is <J Az/(1+z) = 0.029 
and <taz/(i+z) — 0.034, respectively. 

For some fraction, /, of the galaxies, the calculation of the 
photometric redshift fails miserably. These erroneous - or "catas- 
trophic" - redshifts often arise due to confusion between the 
Lyman break (Al = 1215 A) and the Balmer break (Ab = 
4000 A). We use a simple model of catastro phic redshifts base d 
on this confusion and the results presented in lllbert et alJ d2006). 
When the Lyman and Balmer breaks are confused either the con- 
dition 

(1 + £ P )A B = (l + z s )A L , (6) 

or the condition 

(l + z p )A L = (1 + 2 S )A B , (7) 

is fulfilled. If the condition expressed by Eq. (|6]l is fulfilled, 
the photometric redshift is under-estimated. This condition can 
only be fulfilled if z s > 2.3. On the other hand, the photomet- 
ric redshift is over-estimated if Eq. (Q is fulfilled. According 
to Fig. 6c in lllbert et all d2006l) over-estimation of z p is ~ 5 
times more frequent th an under-estimatio n. In our simple model 
a fraction / = 1.9% dllbertet al.ll2006l) of the bright galaxies 
(17.5 < i' AB < 22.5) are assigned catastrophic redshifts. The 
proportion of catastrophic redshifts computed using Eq. (|6]l and 
Eq. dUl is 1:5. For the faint galaxies (22.5 < i' AB < 24) we 
use the same mode l, but the fractio n of catastrophic redshifts is 
higher, / = 5.5% dllbert et alJl2006l) . 

Scatter in the F-J and T-F relations In the conversion be- 
tween galaxy luminosity and halo mass, the Faber-Jackson and 
Tully-Fisher relations are used. The scatter in these relations, 
given by Eq. (01 and Eq. ((5]), are the largest sources of uncer- 
tainty in the calculations. When p cst is estimated this scatter is 
taken into account. We also investigate the effects of system- 
atic errors in the relations between luminosity and mass (see 
Sect. 15.51 ). Missclassification of galaxy type is a related potential 
uncertainty which we have not investigated. However, since the 
missclassification rate for the brightes t galaxies, which are also 
the most important lenses, is likely dDahlenll2008l) to be only a 
few percent, we assume the neglection of this source of uncer- 
tainty to be safe. 

Halo models For the fiducial model, all galaxies are described 
by NFW profiles. When we estimate the uncertainties in the 
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magnification factor we assume incorrect halo models, i.e. SIS 
instead of NFW profiles, for 50% of the lenses. 

4.4. Correlation 

To predict the probability of detecting a correlation between es- 
timated magnifications and Hubble diagram residuals we use 
Monte Carlo simulations. A large number of SNe la data sets are 
simulated and for each data set the linear correlation coefficient, 
r, is computed. The residuals in the Hubble diagram, which is 
built of SNe la corrected for light curve shape and color, de- 
pends on the true magnification factor (computed using the fidu- 
cial model), 

A = -2.51og 10 /z. (8) 

In our simulations we add gaussian noise to the Hubble di- 
agram residuals from intrinsic brightness scatter and measure- 
m ent errors. For the intrinsic dispersion we use the value derived 
bv lAstier et al] j2006). Omt =0.13 mag. To model the measure- 
ment errors in the SNLS we use the following reds hift dependent 
fit to the measurement errors in lAstier et alj d2006h 

J 0.05 mag if z < 0.8 

cr ° IT ~~ \ 0.84z 2 - 1.04^ + 0.34 mag if z > 0.8. { ' 

The total noise in the residuals are given by <7i n t and cr err added 
in quadrature. This noise or other sources of noise could in real- 
ity be non-gaussian. Such non-gaussianity will be a concern for 
techniques that try to detect gravitational lensing of SNe la based 
on th e asymmetry of the distribution of residuals (e.g., Wang 
120051) . However, for this study where we calculate the magnifi- 
cation for each individual supernova, as long as the contributions 
to the residuals are not correlated with the foreground matter, 
the gravitational lensing signal should not be systematically dis- 
torted, only diluted. 

The magnification, in logarithmic units, on the other hand, 
depends on the estimated magnification factor, 

A^ = —2.5 log 10 /West- (10) 

The uncertainties contributing to the noise in the estimated mag- 
nification factors were described in Sect. 14.31 

Since we assume the intrinsic brightness scatter and mea- 
surement errors to be gaussian and these errors dominate, the 
distributions of correlation coefficients are gaussian as well and 
can be described by the mean value, (r), and the standard devia- 
tion, ay. To evaluate the confidence level of the detection of the 
correlation, we compare with what we could expect to measure if 
there is no correlation. For each simulated data set we compute 
the correlation coefficient of the null hypothesis, r nu n> i-e. the 
hypothesis of no lensing. Since both distributions are gaussian, 
the probability of a na detection is given by 

where (r nu n) and oy null refers to the mean and standard deviation 
of the distribution of r nu n, respectively. 

5. Results 

5. 1. Gravitational magnification 

In this section we present some basic properties of the simulated 
gravitational magnification distributions. 



5.1.1. SNLS 

We have simulated the magnification distribution for sources fol- 
lowing the redshif t distribution of the first year SNLS SNe la 
(lAstier et alj|2006l ). The magnification factors used to build up 
this composite distribution were drawn from redshift dependent 
magnification distributions, which were all normalised to fulfil 
the condition (/i) = 1, Although the average magnification fac- 
tor is unity, the average value of A, which depends on the higher 
moments of the magnification distribution, deviates from zero, 
(A) = 0.0011 mag. Of importance is also the median of the 
magnification distribution, A1/2 = 0.0084 mag, which tells us 
that most SNe la in the SNLS data set are slightly demagnified, 
and would consequently appear to be more distant than they re- 
ally are. 

5.1.2. SDSSII 

According to our simulations, the average and the median of 
the composite magnification distribution corresponding to the 
SDSSII SN la survey is (A) = 0.00018 mag and A 1/2 = 
0.00080 mag, respectively. The magnification of most SNe la 
are thus clearly negligible for the SDSSII, which is of course 
a consequence of the shallower redshifts probed by this survey. 
Since the expected effect of gravitational lensing is very small 
for the SDSSIISN, we focus on gravitational lensing of SNLS 
SNe la for the rest of this paper. 

5.2. Shifts in cosmological parameters 

Gravitational lensing can potentially lead to biased cosmological 
results. Here we investigate and quantify these effects for the 
SNLS and show that the effects are indeed rather small. We have 
studied shifts in cosmological parameters for two cases: 

1 . The universe is assumed to be flat and dominated by a cos- 
mological constant. In this case the present matter density, 
fiM is fitted to simulated data sets. 

2. The universe is assumed to be flat and the dark energy equa- 
tion of state parameter, w, is furthermore assumed to be con- 
stant. In this case Qm and w are simultaneously fitted to data 
and constrains from Baryonic Acoustic Oscillations (BAO, 
Eise nstein et al.ll2005h are used. 

A simulated low redshift data set consisting of 44 S Ne la, with 
properties similar to that used by lAstier et al.1 (120061) . is used to 
anchor the Hubble diagram. 

We characterise the shift in S1m and w by the differences, 
Ail M = f^ nscd - n M and Aw = w lcnscd - w , in the best 
fit values obtained with and without gravitational lensing. Our 
results are collected in Table Q] 

Since flux is conserved, effects of gravitational lensing will 
average out for large numbers of SNe la sampling the whole 
magnification distribution. In this study, however, magnitudes 
are used instead of fluxes, as is common in supernova cos- 
mology. When magnitudes are used, the average magnitude is 
shifted compared to the case of no lensing. This shift is caused 
by higher order moments of the magnification factor distribution, 
which are responsible for the fact that (A) 7^ even though 
(/i) = 1. We thus expect a small shift in the cosmological pa- 
rameters solely due to our choice of using magnitudes. For the 
SNLS, where (A) = 0.0011 mag, these shifts are expected to 
be A^ M = -0.0012 for case 1 and Atl M = -0.00030 and 
Aw = -0.0024 for case 2. 
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Table 1. Average and root mean square of distributions of AHm 
and Aw for case 1 and case 2 for different SNLS like data sets 
with and without bright outliers removed. 





Outliers^ 


Case 1 


Case 2 


N 


removed 


(Aft M ) 




(Aw) 




70 


no 


-0.0013 


0.0075 


-0.0020 


0.013 


70 


yes 


-0.0038 


0.010 


-0.0043 


0.013 


500 


no 


-0.0016 


0.0045 


-0.0023 


0.0059 


500 


yes 


-0.0052 


0.0061 


-0.0051 


0.0056 



" Bright SNe la deviating more than 2.5<j from the best fit cosmology. 



Figure Q] shows the shift in £Im which, according to or sim- 
ulations, can be expected for SNLS first year and final data sets. 
The top and bottom panel shows the distribution (shaded his- 
tograms) of A51m for 70 and 500 SNe la, respectively. Here we 
consider the final data set which will be used for cosmology and 
not the smaller subset which could be used to detect lensing. 
For N = 70 the distribution is clearly asymmetric, while for 
N = 500 the distribution is fairly gaussian. The average shift, 
(AQm)> or bias, is very close to the value expected for both 
N = 70 and 7Y = 500 . 

lAstier et"aT] d2006l) estimate Hm obtained from the first year 
SNLS data set to be shifted due to gravitational lensing by at 
most AQm = —0.005. The distribution in Fig.Q]peaks close to 
this value for N = 70. We expect the contribution to the uncer- 
tainty in £Im from gravitational lensing to be of the order 0.008 
or 3% at the la (68.3%) confidence level. Much of the disper- 
sion due to gravitational lensing should, however, already be in- 
cluded in the SNLS analysis via the intrinsic dispersion derived 
from the data themselves. For the final data set, which will con- 
sist of ~ 500 SNe la, the expected contribution to the error in 
SIm from lensing is ^ 2% at the la confidence level. 

Let us now consider the shift for case 2. Figure [2] shows the 
results of our simulations for SNLS. To each simulated data set, 
f^M and w were fitted simultaneously with and without lensing. 
The shaded histograms in the top and bottom panel in Fig. [2] 
show the probability distributions of Aw for N = 500 and N = 
70, respectively. Since the shift in JIm is very small for case 2, 
we only consider Aw. Also for case 2 the expected bias, (Aw), 
is close to what we expect due to the higher order moments in the 
magnification factor distribution (outlined by the vertical lines in 
the figure). 

The most likely shift for N = 70 is Aw « —0.0 05, which is 
half th e maximum shift (Aw = —0.01) estimated bv lAstier et alj 
(2006). The uncertainty in w from gravitational lensing at the lcr 
confidence level for the first year and final SNLS data together 
with BAO constraints is expected to be at the percent and sub- 
percent level, respectively. Most of this uncertainty, as already 
noted, should be included in the intrinsic scatter derived for the 
SNLS. 

5.3. Removal of outliers 

Since the magnification distribution is asymmetric, it is impor- 
tant to sample the high magnification tail. Highly magnified 
SNe la would be outliers in the Hubble diagram and might there- 
fore be removed. 

ISarkar et alj d2008l) recently quantified the shift in the equa- 
tion of state parameter which could be expected due to lensing 
for future large data sets (<; 2000) of SNe la. They also investi- 
gated the effect of removing bright outliers. According to their 
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Figure 1. Shift in £Im due to gravitational lensing. The grey 
shaded distributions show Ail M = ^M nscd - for - 5000 
simulated data sets consisting of 70 (bottom panel) and 500 (top 
panel) SNLS SNe la. The vertical solid line shows the shift 
(AJ7m = —0.0012) expected due to higher order moments of 
the magnification distribution. The hatched histograms show the 
shift which would be the result of removing bright 2.5cr outliers 
from the Hubble diagram. 



simulations, the shift in w increases from ~ 0.5% to ~ 0.8% 
when bright SNe la deviating by more than 2.5<r are removed. 

We have performed a simple exercise to investigate the ef- 
fects of removing highly magnified SNe la from the Hubble dia- 
gram. The hatched histograms in Figs.[T]and|2]show the effect of 
removing bright 2.5er outliers. Approximately 1% of the SNe la 
are removed by this cut. The results of the exercise are given in 
Table □ 

Let us first discuss the effect of removing bright outliers for 
case 1 . If bright outliers are removed, the peak of the distribution 
of A17m is shifted to a lower value. Moreover, the distribution 
broadens and become more gaussian. Since the high magnifica- 
tion SNe la have been removed, most SNe la appear to be dim- 
mer and consequently the data favour a smaller value of J7m- The 
removal of the asymmetric high magnification tail also explains 
why the distribution of AHm becomes more gaussian. The av- 
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Figure 2. Shift in w due to gravitational lensing. The grey 
shaded distributions show Aw — w — ^ lensed for ~ 5000 sim- 
ulated data sets consisting of 70 (bottom panel) and 500 (top 
panel) SNLS SNe la. The vertical solid line shows the shift 
(Aw = —0.0024) expected due to higher order moments of 
the magnification distribution. The hatched histograms show the 
shift which would be the result of removing bright 2.5cr outliers 
from the Hubble diagram. Constraints from BAO were used in 
the cosmology fits. 



erage shift in fl^i increases with approximately a factor 3 when 
the outliers are removed. 

The peak of the distributions shift towards lower values also 
for case 2. The shift in w increases by a factor approximately 2 
when the bright outliers are removed. However, the width of the 
distributions hardly change at all. 

5.4. Detecting a magnification correlation 

Since amplified and de-amplified SNe la should be brighter and 
fainter than average, respectively, we expect a correlation be- 
tween the residuals in the Hubble diagram and the estimated 
lensing magnifications, A^. Clearly, intrinsic brightness scatter 
among the SNe la as well as measurement errors add noise to the 
gravitational lensing signal. The ability to detect the correlation 
therefore depends on the quality and size of the data set. 

Our simulations demonstrate that the quality of the SNLS 
data is sufficient for a detection of a correlation at high confi- 
dence. Figure [3] shows the probability of detecting the correla- 
tion at lcr (68.3%, solid curve), 2a (95.4%, dashed curve), and 
3a (99.7%, dotted curve) confidence level as a function of the 
number of SNe la. For these simulations we have used the dis- 
tribut ions of redshifts and errors from the first year of SNLS 
data ( lAstier et~a l. 2006). If the magnification factor could be es- 
timated for 450 SNe la, the probability to detect the correlation 
at the 3cr confidence level is 97.6%. A firm detection of the cor- 
relation between magnifications and SN la brightnesses is thus 
very likely to be found from the final SNLS data. If the magni- 
fication factor could be predicted exactly, the probability would 
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Figure 3. Probability to detect a correlation between estimated 
magnification and SN la brightness at different confidence lev- 
els, P(no), as a function of the number, N, of SNe la. The solid, 
dashed, and dotted curves shows the probability to detect a cor- 
relation at the la, 2<7, and 3cr confidence level, respectively. The 
simulated data sets have the same distributions of redshifts and 
errors as the first year SNLS data. 



instead be 99.9%. Even in this best case scenario uncertainties 
in the Hubble diagram residuals would affect our possibility to 
find the correlation. 

Increasing the intrinsic brightness dispersion from cr; n t = 
0.13 mag to <7i, lt = 0.15 mag and <j mt = 0.17 mag would 
decrease the probability to make a 3cr detection of the lensing 
signal from 97.6% to 96.1% and 90.2%, respectively. The pos- 
sibility to detect a lensing signal at high confidence level is thus 
not critically dependent on the intrinsic dispersion. 

Our chances to make a high confidence level detection of the 
correlation using the final SNLS data looks promising even if our 
estimates of the lensing masses would be wrong. The probability 
to make a 3<7 detection is J> 95% even when the lensing masses 
are over-estimated or under-estimated by ~ 50%. For the final 
SDSSIISN data set consisting of ~ 300 SNe la, on the other 
hand, we predict the probability of a 3cr detection to be only 4%. 



5.5. Constraints on halo masses 

The gravitational lensing signal provided by the correlation be- 
tween magnifications and SN la brightnesses should also give 
some constraints on masses of the haloes of the lensing galax- 
ies. We have performed simulations where the mass of each halo 
has been multiplied by a factor F. The increase or decrease in 
halo masses is compensated for by changes in the smoothness 
parameter to ensure that the total mass in the Universe sums up 
to the global value. To each simulated data set, a linear relation, 
A = A + BA^, was fitted. Figure [4] shows how the slope, de- 
scribed by the parameter B, changes with the value of F. The 
dashed horizontal line corresponds to B = 1, which is the value 
which would be recovered if the magnification factors were ex- 
actly known. If obvious outliers in the magnification-residual di- 
agram are removed, the distribution of the B parameter is gaus- 
sian. The inner and outer error bars in Fig. [4] correspond to lcr 
and 2cr errors. The number of SNe la was N = 450, however 
the standard deviation scales roughly as 1/yN, 
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Figure 4. Slope, B, in the magnification-residual diagram as a 
function of the factor F which all halo masses have been multi- 
plied by. The vertical dashed line shows the value expected for a 
correlation, i.e. 5 = 1. The inner and outer errorbars correspond 
to la and 2a confidence level, respectively. 
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According to Fig. [4] the slope in the magnification-residual 
diagram varies with F, which should allow us to probe the nor- 
malisation of the halo masses. From the figure we conclude that 
the normalisation of the halo masses could probably be con- 
strained with the final SNLS data set at the la and 2a confidence 
level with ~ 30% and ~ 60% accuracy, respectively. 



6. Conclusions 

We have simulated the effect of gravitational lensing in two ma- 
jor ongoing supernova surveys. For the relatively nearby SDSSII 
supernova search, the effect of gravitational lensing will be 
small. For the more distant supernovae in the SNLS survey, we 
predict that the signal from gravitational lensing will be observed 
with high confidence. Our simulations indicate that a correlation 
between Hubble diagram residuals and magnification for indi- 
vidual supernovae will be present at high (at least 3tr) signif- 
icance level. This could be used both to somewhat reduce the 
scatter in the Hubble diagram, and to learn about the proper- 
ties of the lensing material. A project to investigate this effect in 
the SNLS data is underway. We also note that the prospects of 
using weak lensing of supernovae to constrain the matter distri- 
bution in the Universe with future surveys such as SNAP look 
very promising. A satellite like SNAP would provide not only 
high redshift SNe la, but also deep observations in many filters 
which would allow reliable photometric redshifts of the galaxies 
to be obtained. 
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